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ABSTRACT 

In the near future, observations of the cosmic microwave background (CMB) 
anisotropics will provide accurate determinations of many fundamental cosmological 
parameters. In this paper, we analyse degeneracies among cosmological parameters to 
illustrate some of the limitations inherent in CMB parameter estimation. For simplic- 
ity, throughout our analysis we assume a cold dark matter universe with power-law 
adiabatic scalar and tensor fluctuation spectra. We show that most of the variance 
in cosmological parameter estimates is contributed by a small number (two or three) 
principal components. An exact likelihood analysis shows that the usual Fisher ma- 
trix approach can significantly overestimate the errors on cosmological parameters. 
We show that small correlated errors in estimates of the CMB power spectrum at 
levels well below the cosmic variance limits, (caused, for example, by Galactic fore- 
grounds or scanning errors) can lead to significant biases in cosmological parameters. 
Estimates of cosmological parameters can be improved very significantly by applying 
theoretical restrictions to the tensor component and external constraints derived from 
more conventional astronomical observations such as measurements of he Hubble con- 
stant, Type la supernovae distances and observations of galaxy clustering and peculiar 
velocities. 



1 INTRODUCTION 



A number of investigations have shown that high precision 
measurements of the cosmic microwave background (CMB) 
anisotropics can be used to determine many cosmological pa- 
rameters to unprecedented precision (Jungman et al. 1996; 
Bond, Efstathiou and Tegmark 1997, hereafter BET97; Zal- 
darriaga, Spergel and Seljak 1997, hereafter ZSS97). These 
parameters include the amplitudes and spectral indices of 
the scalar and tensor fluctuations predicted by inflationary 
models (e.g. Knox and Turner 1994, Knox 1995, Lidsey et al. 
1997, Souradeep et al. 1998, Copeland et al. 1998), the den- 
sities of the various components of the mass density of the 
Universe, and constraints on the geometry of the Universe. 

The prospect of measuring such cosmological parame- 
ters to high precision has formed an important part of the 
scientific case for two approved satellite missions, NASA's 
MAP satellite (Bennett et al. 1996) and ESA's Planck Sur- 
veyorQ(Bersanelli et al. 1996). In addition, ground based and 
balloon experiments are already beginning to set useful con- 
straints on some of these parameters (e.g. Bond and Jaffe 
1997, Hancock et al. 1997, Lineweaver and Barbosa 1998a,b, 
Webster et al. 1998). 

* Formerly known as COBRAS/SAMBA 



In this paper, we investigate degeneracies among cos- 
mological parameters, i.e. to what extent can we disen- 
tangle one parameter from another using measurements of 
the CMB anisotropies alone (see also Bond et al. 1994). 
We assume throughout that the primordial fluctuations are 
adiabatic and Gaussian, as expected in most inflationary 
models (e.g. Linde 1990). These assumptions are physically 
well motivated and allow accurate calculations of the CMB 
anisotropies as a function of cosmological parameters via 
fast numerical solutions of the relativistic Boltzmann equa- 
tion. The accuracy of computations of CMB anisotropies in 
topological defect models has greatly improved recently (see 
Turok, Pen and Seljak, 1997, Albrecht, Battye and Robinson 
1997) and tend to disfavour the simplest versions of defect 
models. However, these calculations cannot yet be done to 
sufficient precision for the type of analysis presented in this 
paper. Nevertheless, some qualitative features of the dis- 
cussion presented here are likely to be applicable to defect 
models and to other classes of models such as isocurvature 
theories. 

We distinguish between a nearly exact degeneracy 
(which we call the 'geometrical' degeneracy) and other pa- 
rameter degeneracies. The geometrical degeneracy (BET97, 
ZSS97) leads to near identical CMB anisotropies in uni- 
verses with different background geometries but identical 
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matter content. In linear perturbation theory, observations 
of the primary CMB anisotropies cannot break the geomet- 
rical degeneracy, no matter how precise the experimental 
measurements. The geometrical degeneracy thus imposes 
fundamental limits on measurements of the curvature of 
the Universe and the Hubble constant derived from CMB 
anisotropy measurements. Of the remaining parameter de- 
generacies, some are extremely sensitive to the accuracy of 
the CMB anisotropy measurements whilst others are not. In 
this paper, we give a physical description of the causes of 
some of the parameter degeneracies and of the fundamental 
limits imposed by CMB anisotropy measurements. 

This paper is laid out as follows. In Section ^| we sum- 
marise the formalism based on the Fisher matrix (Jungman 
et al. 1996, Tegmark, Taylor and Heavens 1997) which is 
assumed in most discussions of parameter estimation from 
the CMB. Throughout this paper, we use approximate ex- 
perimental parameters for the MAP and Planck satellites 
to provide numerical examples of what might be achieved in 
the next decade or so. The experimental parameters that we 
use are summarised in Section ^. In Section ^, we discuss the 
geometrical degeneracy and we show that with an appropri- 
ate change of variables, the Fisher matrix can be applied 
to give accurate results for cosmological models of arbitrary 
curvature. In Section ^, we specialise to spatially flat cosmo- 
logical models and analyse degeneracies among other cos- 
mological parameters. We show that a principal component 
analysis provides a powerful technique for analysing degen- 
eracies among parameters. For experiments such as MAP 
and Planck, that sample the CMB anisotropy power spec- 
trum beyond the first Doppler peak, the first few principal 
components dominate the expected variances of most cos- 
mological parameters. We test the accuracy of the Fisher 
matrix approach by performing an exact likelihood analysis 
in a subspace defined by the two most poorly determined 
principal components. Some of the degeneracies among cos- 
mological parameters can be understood in terms of the lo- 
cations of the Doppler peaks and of the height of the first 
Doppler peak. In Section [|, we compare the accuracies of pa- 
rameter estimates from the improved version of MAP with 
those from Planck. We distinguish between 'strength' and 
'statistical significance' in the analysis of CMB power spec- 
tra and discuss the effects of possible systematic errors that 
might arise, for example, through inaccurate subtraction of 
point sources, Galactic foregrounds and map reconstruction 
errors. 

Section [| discusses briefly how some of the degeneracies 
among cosmological parameters can be removed by using 
other types of observation, e.g. measurements of the Hub- 
ble constant, the age of the Universe and constraints on 
the luminosity distance from observations of distant Type 
la supernovae. Our conclusions are summarised in Section 
|(| Throughout this paper we use CMB power spectra com- 
puted with the CMBFAST code developed by Seljak and 
Zaldarriaga (1996). 



FISHER MATRIX AND EXPERIMENTAL 
PARAMETERS 



on the celestial sphere can be decomposed in a spherical 
harmonic expansion, 
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where each coefficient a^ m is statistically independent and 
Gaussian distributed with zero mean and variance 



Ci = (|o/, TO | ) • 



(2) 



The power spectrum, Ci, completely specifies the statisti- 
cal properties of the primary CMB temperature fluctuations 
if the initial fluctuations are Gaussian. The anisotropy of 
Thomson scattering results in a small net linear polariza- 
tion of the CMB anisotropies (Kaiser 1983, Bond and Ef- 
stathiou 1984). For near scale- invariant spectra, this linear 
polarization amounts to a few percent of the total temper- 
ature anisotropy and should be measurable in future CMB 
experiments such as MAP and Planck. Polarization intro- 
duces three additional power-spectra to fully characterise 
the temperature anisotropies. These define the 'electric' and 
'magnetic' components of the polarization pattern and the 
cross correlation of the polarization pattern with the temper- 
ature anisotropies (see e.g. Seljak and Zaldarriaga 1996, Hu 
and White 1997). We do not analyse polarization in any de- 
tail in this paper. The effects of polarization measurements 
on cosmological parameters has been discussed recently by 
ZSS97. Their results show that polarization measurements 
can improve the accuracy of most cosmological parameters 
by relatively modest factors, with two important exceptions, 
distinguishing between scalar and tensor modes and deter- 
mining the redshift at which the intergalactic medium was 
reionized. The role of polarization will be discussed further 
in Section ^| 

We follow the notation of BET97 and assume that a 
cosmological model is specified by a set of parameters s . 
Let P(s\prior) be the prior probability distribution of the 
parameters s and C(s) the likelihood function defined by the 
experiment. The probability distribution of the parameters 
s, taking into account the new experimental information, is 
given by Bayes theorem P(s) oc £.(s)P(s\prior). If the errors 
As = s — so about the target model parameters so are small, 
an expansion of ln£ to quadratic order about the maximum 
leads to the expression, 



£ w C m exp 



2 X/ Fi 3 SSiSs 3 
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(3) 



where F,j is the Fisher matrix, given by derivatives of the 
CMB power spectrum with respect to the parameters s 
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(4) 



The quantity ACi in equation ^| is the standard error 
on the estimate of Ct. For an experiment with N fre- 
quency channels (denoted by subscript c), angular resolu- 
tion (9 C ) f-uihm and sensitivity (a c ) p i x per resolution element 
((9c)fwhm x (6 c ) fwhm pixel) sampling a fraction f sky , 



We assume that the initial fluctuations are adiabatic and 
Gaussian. The temperature anisotropies AT /T of the CMB 
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w = w c , B e = Bciw c /w, (5b) 

c c 

w c = (a c ,pixVc, P ix) , o ct ^ e j (5c) 

(see Knox 1995 and BET97), where we have assumed that 
the experimental beam profile £> c is Gaussian with width 
l s = \/81n2(# c ) J^ hm . If we assume a uniform prior, as we will 
do throughout this paper, then the covariance matrix M = 
(<5s5s^) is the inverse of the Fisher matrix F. The standard 
deviation of a parameter Si, marginalized over uncertainties 

1/2 

in the other parameters, is given by Oi = M M . 

Table || lists the resolutions, @fwhm> sensitivities Opi X 
and the noise power parameters w~ for the various fre- 
quency channels that we have adopted. The parameter set 
labelled OMAP ('original MAP') is for the original specifica- 
tions of the MAP satellite, as given in the proposal of Bennet 
et al. Parameters for the current MAP design are listed un- 
der the heading CMAP ('current MAP'). These include a 
significant improvement in the angular resolution at all fre- 
quencies compared to the original design specifications. We 
have listed the parameters for the 40, 60 and 90 GHz chan- 
nels; the angular resolutions of the lower frequency MAP 
channels are so much poorer (0.93° at 22GHz and 0.68° at 
30GHz) that they carry very little weight in the estimation of 
cosmological parameters and so are ignored here. The main 
purpose of these lower frequency channels is to monitor the 
free-free and synchrotron contributions from the Galaxy. 

For Planck, we have adopted parameters for the four 
lowest frequency channels for the high-frequency instrument 
(HFI), as in the current design submitted to ESA (Puget et 
al. 1998). The Planck payload is in a state of active devel- 
opment and may change as the design of the Planck instru- 
ments is refined. However, the parameters listed in Table ^ 
give an accurate indication of what Planck is intended to 
achieve. The major difference in the current Planck payload 
compared to that described in the Phase A report is the 
addition of polarization sensitivity in the HFI instrument at 
143, 217 and 545 GHz and a much improved performance 
of the low-frequency instrument (LFI) , which covers the fre- 
quency range 31 - 100 GHz. The proposal to cool the LFI to 
20K is expected to improve its sensitivity to within a factor 
of two (or better) of that of the HFI. The parameters listed 
in Table || for Planck are thus likely to underestimate the 
total sensitivity achievable by summing over all frequency 
channels. 

Evidently, there are some uncertainties in the experi- 
mental parameters, particularly for Planck. We emphasise, 
therefore, that the numbers given in this paper are meant 
to be illustrative of the degeneracies among cosmological 
parameters resulting from experiments with the character- 
istics of MAP and Planck, rather than a precise analysis of 
these satellites. A detailed analysis should include system- 
atic errors arising from Galactic and extra-Galactic fore- 
grounds, scanning strategy and map reconstruction, and 
power-spectrum estimation. A considerable amount of work 
along these lines has been done already {e.g. Bouchet, Gis- 
pert and Puget 1995, Tegmark and Efstathiou 1996, Bouchet 
et al. 1997, Hobson et al. 1998, Delabrouille, 1998) and will 
be developed in the future using more precise models of the 
satellites. In this paper, we assume that the CMB power 
spectrum can be measured to the theoretical variance given 
by equation ^3j ) over half of the sky, f s k y = 0.5. We there- 



fore ignore any non-primordial contributions to the CMB 
anisotropies (except in Section |H|) , assuming that these can 
be subtracted to high accuracy using the frequency informa- 
tion available with both satellites. 

The three sets of parameters in Table ^ are also useful 
for presentational purposes, since they help illustrate some 
of the key points of this paper. The OMAP specifications 
are close to a threshold at which there are large degenera- 
cies between estimates of cosmological parameters. These 
specifications have been retained for their pedagogical value, 
even though the actual MAP performance is expected to be 
much better. As we show in Sections |^ and |i|, the CMAP 
parameters (primarily the higher angular resolution) lead 
to a significant improvement in cosmological parameter es- 
timates. Of course, this is the main motivation for improv- 
ing the angular resolution of MAP compared to the origi- 
nal design specification. The parameters for Planck define 
an experiment that is almost 'cosmic- variance limited', i.e. 
where random experimental errors are less important than 
the sampling errors arising from the fact that we can ob- 
serve only one realisation of microwave sky. In summary, we 
will see in later sections that degeneracies among cosmolog- 
ical parameters are severe for OMAP, the results for Planck 
are close to the best that can be achieved from observations 
of the primary CMB anisotropies alone, and the results for 
CMAP are intermediate between these two cases. 



3 GEOMETRICAL DEGENERACY 
3.1 Cosmological parameters 

We specify a spatially flat target cold dark matter (CDM) 
model defined by the following parameters: 

• Spectral indices n s and nt for the scalar and tensor com- 
ponents. The target model has n s = 1 and n t = 0. 

• The overall amplitude of the CMB anisotropy spectrum 
defined by the average band power {£(£ + l)Ce/(2n)) 1 ^ 2 over 
the range of multipoles £ m ax accessible to the experiment. 
We denote this amplitude, relative to a COBE normalized 
model, by Q. The target model has Q = 1. 

• The ratio r of the tensor and scalar components r = 
Cl/Cf. The target model has r = 0.2.f] 

• The cosmological densities, Qi , of the various constituents 
of the universe, relative to the critical density of the Ein- 
stein -de Sitter model. We include the contribution from 
baryons fib, cold dark matter £7 C , a cosmological constant 
Qa = A/(3Hq) (where Ho is Hubble's constant), and a den- 
sity parameter that characterises the spatial curvature of the 
Universe: 

f2fe = 1 — Of, — Q c — Qa ■ (6) 

The CMB anisotropies are determined directly by physical 
densities, uji, which are related to the density parameters Qi 

t In inflationary models with an exactly uniform rate of acceler- 
ation the ratio r and the tensor and scalar spectral indices are 
related via r — 7nt = 7(1 — n s ) (see e.g. Davis et al. 1992, 
Liddle and Lyth 1992, Bond et al. 1994). The parameters of the 
target model are slightly inconsistent with these relations, though 
this is of no consequence for any of the subsequent results in this 
paper. 
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Table 1: Experimental Parameters 



OMAP CMAP PLANCK 



v (GHz) 


40 


60 


90 


40 


60 


90 


100 


150 


220 


350 


@ fwhm 


32' 


2.3' 


17.4' 


28' 


21' 


12.6' 


10.7' 


8' 


5.5' 


5' 


<xp»/icr 6 


7.6 


10.8 


16.3 


6.6 


12.1 


25.5 


1.7 


2.0 


4.3 


14.4 




4.9 


5.4 


6.8 


2.9 


5.4 


6.8 


0.028 


0.022 


0.047 


0.44 



by LOi = Qih 2 , where h is Hubble's constant Ho in units of 
100 kms - Mpc - . We therefore use the physical densities U)i 
to specify the target model. The Hubble constant then be- 
comes a secondary parameter that is fixed by the constraint 
equation 



■■ UJb + LO c + LU k + WA • 



(7) 



The target model has u>b = 0.0125, uj c = 0.2375, ujk = 0, 
luk = 0, i.e. Q, b = 0.05, D.c = 0.95, fl k = 0, fi A = and h = 
0.5. The advantages of using the densities cui for parameter 
estimation have been discussed previously by BET97, and 
will become obvious in the next two sections. 

The target model is therefore specified by 8 parameters, 
4 parameters defining the spectrum of fluctuations (two am- 
plitudes and two spectral indices) ; and 4 density parameters 
specifying the matter content and geometry of the universe. 
This is a relatively small parameter set, but has been chosen 
to illustrate some key points concerning parameter degen- 
eracies. There would be little value in including, for example, 
additional matter components such as hot dark matter, or 
relativistic matter, since this would complicate the discus- 
sion of degeneracies without adding anything substantially 
new. 

We have not included any parameters to specify reion- 
ization of the intergalactic medium, e.g. an optical depth 
to Thomson scattering, r op t. This is because in the absence 
of polarization measurements, the effects of reionization are 
almost completely degenerate with a change in the ampli- 
tude of the fluctuations (i.e. the product Qe Topt is equal 
to a constant). With polarization information, it is possible 
to break the degeneracy between Q and T opt , provided r op t 
is large enough (see ZSS97). The measurement of polariza- 
tion also offers a possibility of distinguishing between tensor 
and scalar modes (Crittenden and Turok 1995, Seljak 1997a, 
Kamionkowski and Kosowsky 1998). If the polarization pat- 
tern can be reliably decomposed into electric and magnetic 
components, it is possible to determine the relative ampli- 
tudes of the scalar and tensor modes to much higher accu- 
racy than is possible from observations of the total CMB 
anisotropy alone. Such a measurement would have a large 
qualitative impact on the parameter estimates described in 
the next two sections, significantly reducing the degenera- 
cies between the tensor parameters and other cosmological 
parameters. However, the tensor contribution to the polar- 
ization signal is very small and concentrated to low I values, 
so may be difficult to measure against contaminating polar- 
ization signals of Galactic origin. 



3.2 The angular diameter distance degeneracy 

Two physical scales play a crucial role in determining the 
shape of the CMB anisotropy spectrum. These are the sound 
horizon r s (a r ) at the time of recombination, and the angu- 
lar diameter distance cLa to the last scattering surface. (The 
quantity a r is the cosmological scale factor at recombination, 
normalized to unity at the present epoch). The sound hori- 
zon defines the physical scales for the Doppler peak struc- 
ture that depends on the physical matter densities u>b and 
uj c , but not on the value of the cosmological constant or 
spatial curvature since these are dynamically negligible at 
the time of recombination. The angular diameter distance 
depends on the net matter density, uj m = u>b + uj c , u)k and 
uja and maps the physical positions of the Doppler peaks 
to peaks in the angular power spectrum Ci as a function 
of multipole I. This mapping leads to a nearly exact 'geo- 
metrical' degeneracy among cosmological parameters, as has 
been noted by BET97 and ZSS97. Models will have almost 
exactly the same CMB anisotropy spectra at high multipoles 
if they satisfy the following criteria: 
[1] they have identical matter densities u)b and ui c ; 
[2] identical primordial fluctuation spectra; 
[3] identical values of the parameter^ 



K = 



772 



sinh 



,1/2, 



where 



V = 



[uj m a + uj k a? + oj A a 4 + Lu Q a 1 - 3u, Q] 1 / 2 ' 



(8a) 



(8b) 



Here ujq refers to a component of the energy density that 
obeys an equation of state with fixed pressure to density 
ratio, wq — p/p, as might arise with scalar fields which 
are important at late times in the universe, e.g., Ratra 
and Peebles 1988 and Caldwell, Dave and Steinhardt 1998 
(who have dubbed such phenomena quintessence). Note that 
Wq = — 1/3 gives the same contribution as ui k gives in y, but 
the geometrical aspects of the u>k angle-distance relation are 
absent, yielding a behaviour more similar to that of uia- 
Since it therefore does not add qualitatively new features 



t Throughout this paper we restrict our analysis to models with 
negative spatial curvature (positive Q k )- There is no fundamental 
reason for imposing this restriction, since it is straightforward to 
generalize expressions such as equation (^) to models with posi- 
tive curvature. However, all of the numerical examples described 
in this paper have been computed with an updated version of the 
CMB anisotropy code of Seljak and Zaldarriaga (1996) which is 
restricted to models with negative spatial curvature. 
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0.1 0.2 0.3 0.4 0.5 



fi k 0) k 

Figure 1. Figure |l]a shows degeneracy lines of constant 1Z in the QA _ ^fc plane. The value of TZ is given next to each line. In computing 
this figure, the matter density parameter is fixed by the constraint f2 m = 1 — £\- — ^A- Figure hp shows lines of constant 1Z in the 
ujfr-Wk plane for a universe with u) m = 0.250. The five dots in each figure show the loc atio ns of the five models with nearly degenerate 
Cf spectra plotted in Figure 0. The target model with the parameters given in Section 3.1 is located at the origin in each panel. 



to the discussion we ignore uiq terms in what follows. Equa- 
tions ( |8aj|8b| ) are special cases of more general expressions for 
the locations of the Doppler peaks given in Section ^. Pro- 
vided the above conditions are satisfied, the only differences 
between the CMB power spectra of models with different 
values of Qa and fife will be those at low multipoles £ ^ 30 
arising from the late-time Sachs-Wolfe effect (see Bond 1996, 
Hu, Sugiyama and Silk 1997, and references therein). 

This is illustrated in Figures jj] and ^. Figure 1 shows 
lines of constant 1Z in the fiA - fife plane (Figure |l]a) and in 
the uja - iOk plane (Figure |l|b), for a fixed value of u) m = 0.25. 
Figure ^a is the more useful, because the degeneracy lines 
are fixed by the two parameters fiA and fife , with the matter 
density specified by the constraint fi m = 1 — fife — fiA- The 
degeneracy lines in the space defined by the physical densi- 
ties u! m , ujk and uja require specification of all three densities 
or, alternatively, two densities and a value of the Hubble con- 
stant (see equation |^) . We have therefore assumed the value 
of u> m of our target model in computing Figure |l|b. 

The five dots in each panel of Figure [l] show the loca- 
tions of models with nearly degenerate Ce spectra. These 
models have the same matter densities as the target model, 
u>b — 0.0125 and u> c = 0.2375, and identical scalar spectral 
index n s — 1 (we compute only the scalar component for this 
comparison) . The five models have different values of uik and 
uja chosen to lie along the degeneracy line of 1Z = 1.94 of the 
target model, as shown in the figures. The resulting CMB 
power spectra are plotted in Figure |5| 

As expected, the power spectra plotted in Figure ^ are 
indistinguishable except at small multipoles, I ,$ 30. The 
lower panel compares the residuals ACe / Ce for these models 
compared to the spatially flat target model. The two unla- 
belled dashed lines show the error on the residuals arising 
from cosmic variance alone, ACe/Ce = ±{2/{2£ + 1)) 1//2 . 
There are several points worth noting from this figure: 

[1] The differences between the models at high multipoles 
(i ^ 200) are small and arise from numerical errors. These 



include small errors in the computation of the degeneracy 
parameter 1Z, but arise predominantly from numerical errors 
in the CMBFAST code. 

[2] Although the numerical errors at high multipoles are 
much smaller than the cosmic variance errors, they accu- 
mulate in the computation of the Fisher matrix (equation ^ 
and can erroneously break the geometrical degeneracy. Thus, 
numerical errors in what should be an exact degeneracy can 
lead to extremely misleading results, in which the errors in 
the density parameters uja and lo^ appear weakly correlated 
and extremely small. This problem is particularly acute for a 
Planck-type experiment in which the Ce estimates are nearly 
cosmic variance limited to high multipoles. 

[3] The differences in the power spectra at low multipoles, 
I ^ 30, arise from a real physical effect - the late time Sachs- 
Wolfe effect - and are not caused by numerical errors. How- 
ever, the differences are smaller than the cosmic variance 
limits and so it is impossible to differentiate between these 
models at a high significance level from measurements of the 
CMB anisotropies alone. (This point is quantified in detail 
in Figure ^| below). 

[4] The power spectra for the low density cosmological mod- 
els assume a particular model of the origin of the primordial 
fluctuation spectrum (in this case, power-law fluctuations 
in the potentials <E> and ^). If we were to allow complete 
freedom to adjust the shape of the power spectrum on large 
spatial scales, then we could compensate for the differences 
caused by the late-time Sachs Wolfe effect. The existence of 
the late-time Sachs- Wolfe effect cannot be used to distin- 
guish between different cosmological models unless specific 
assumptions are made concerning the shape of the fluctua- 
tion spectrum. 

The futility of breaking the geometrical degeneracy 
from linear CMB anisotropies is illustrated in Figure ^[ In 
this figure, we have computed a grid of scalar Ce spectra 
varying the parameters uja and uik , whilst fixing the remain- 
ing parameters, uit, u> c and n s to those of the spatially flat 



© 0000 RAS, MNRAS 000, 000-000 



6 G. Efstathiou and J.R. Bond 




Figure 2. The upper panel shows the scalar power spectrum for the nearly degenerate models shown by the dots in Figures [lj. The 
target model has loa = = 0, u>\, = 0.0125 and uj c = 0.2375 and n s = 1. The lower panel shows the residuals of these open universe 
models with respect to the spatially flat target model. The numbers next to each curve give the value of for each model. The two 
lines with long dashes show the standard deviation of the residuals from cosmic variance alone. 



target model. Figure |3j shows the exact likelihood ratio, 
C \ (C e (s) - C e (s )) 2 



21n 



E 



(AC, 



(9) 



for an experiment with the OMAP parameters listed in Ta- 
ble [| and £ ma x set to 1000. The contours show approximate 
1, 2 and 3a contours in the uja and ujk plane assuming all 
other parameters are known (i.e. we show contours where 
— 2\n(£./Cmax) equals 2.3, 6.2 and 11.8 corresponding to 
68%, 95% and 99.7% confidence regions if -2\sx{C/ C ma x) 
is x 2 distributed with two degrees of freedom) . The la con- 
tours are broken up artificially by the discrete size of the grid 
over which the Ce's were computed numerically. The 2a con- 
tour limits uja to be less than about 0.7 and uik to be less 
that about 0.1. The contours follow the line of constant 1Z 
computed in Figure [l]b and show that measurements of the 
CMB anisotropics alone cannot remove the geometrical de- 
generacy between uja and uik ■ Improving the accuracy of the 
CMB measurements simply causes the likelihood contours to 
narrow around the degeneracy lines. Although this has im- 
portant consequences for the estimation of other cosmolog- 
ical parameters, improving the accuracy of the experiment 
does not help to remove the geometrical degeneracy. For 
all practical purposes, the geometrical degeneracy in linear 
theory is exact and can be removed only by invoking other 
constraints on the geometry of the Universe {e.g. as derived 
from Type la supernovae, Perlmutter et al. 1997, 1998, see 



Section ^), or from the effects of gravitational lensing on 
the CMB (Stompor and Efstathiou 1998, Metcalf and Silk 
1998). 



3.3 Parameter estimation including spatial 
curvature 

In this Section, we describe how the Fisher matrix formalism 
can be adapted to handle the geometrical degeneracy. This 
can be achieved by redefining the variables s so that one of 
the derivatives dCi/dsi is computed along the direction of 
the degeneracy (i.e. along lines of constant 1Z, as in the ex- 
amples plotted in Figure ^). The derivative of Ce along the 
degeneracy line can then be set to zero at high multipoles 
by construction, eliminating the effects of numerical errors 
in computing the elements of the Fisher matrix. This tech- 
nique has been used previously by us in BET97, and gives 
extremely stable results. 

To illustrate the method, we analyse the errors in the 
cosmological parameters for sm all v ariations around the tar- 

The target model has 
and so the derivatives of 



get model defined in Section 3.1 
Qk = 0, S7a = and fi m = 1 
1Z are 

/ &R_\ _ I &]Z\ _ _1 



(10a) 



where the subscript t denotes that the derivatives are com- 
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Table 2: la errors on estimates of 
cosmological parameters 



0.15 



Figure 3. Likelihood ratio contours in the luj^-uj^ plane for mod- 
els containing only scalar modes. The models have fixed values 
of uif, = 0.0125 and oj c = 0.2375. The likelihood ratio for each 
model has been computed assuming the experimental parameters 
of OMAP (Table ^) and our standard spatially flat target model. 
The contours show where — 2\n(C/ C m ax) has values of 2.3, 6.2 
and 11.8 corresponding approximately to ltr, 2tr and 3c contours 
assuming all other parameters are known. The full line shows the 
degeneracy curve with the same value of TZ as the target model 
(cf Figure hp) and the dashed line shows the T aylo r series ap- 
proximation to this curve, ^ = lui^ (see Section |3- 3|) . 



puted at the parameter values of the target model. For small 
variations around the target model, the degeneracy lines in 
TZ follow 



S£l k = -SQ.A, 



5tok = -Suja- 



It is therefore useful to define an auxiliary density 
7uJk 



UJD 



UJA 



(10b) 



(11) 



to replace uj k in the Fisher matrix analysis (where the sub- 
script, D, denotes 'degeneracy') . For small variations around 
the target model, the constraint ujd — defines the degener- 
acy line TZ = constant. The CMB spectra plotted in Figure ^ 
all have ujd ~ and identical values of ui b and ui c . They can 
therefore be used to define the derivative dCe/duJA, which 
can be set to zero at, for example, £ ^ 200 to suppress the 
effects of numerical errors. The results of such an analysis for 
the experiments summarized in Table ^ are listed in Table 



The results from Table 3.3 show that the la errors in 



uja. are close to unity and extremely insensitive to the ex- 
perimental details. This is in agreement with the likelihood 
analysis of Figure |l| The large error in lua thus reflects 
the geometrical degeneracy. In contrast, the error on the 
auxiliary density parameter ujd decreases from 0.17 for the 
OMAP experimental parameters, to 0.014 for the Planck 
parameters. The error in ujd is set by the accuracy to which 
the parameter TZ can be determined by the experiment, i.e. 
by the ability of the experiment to fix the positions of the 
Doppler peaks. An experiment such as Planck samples the 



Parameter 


OMAP 


CMAP 


PLANCK 


5u b /ui b 


0.075 


0.042 


0.0064 


5u) c /u) c 


0.24 


0.13 


0.020 


5lo d 


0.17 


0.097 


0.014 


5ui A 


1.23 


1.20 


1.09 


SQ 


0.013 


0.0041 


0.0013 


8r 


0.41 


0.27 


0.15 


8n s 


0.059 


0.034 


0.0065 


5n t 


1.015 


0.998 


0.911 




0.2 0.4 0.6 0.8 1 

n 



Figure 4. Results of Monte-Carlo simulations of parameter es- 
timation as described in the text assuming the parameters of 
OMAP as given in Table ^ The panels show correlations be- 
tween various pairs of parameters. The lines in the lua-uj i. an d 
f!A-f2fc plane show the degeneracy lines given by equation ( jlObj ). 
The line in the h-QA plane shows the degeneracy line given by 
equation (113) . 



entire Doppler peak structure of Ci and so can determine TZ 
(and hence ujd) to extremely high accuracy. 

Another important consequence of the geometrical de- 
generacy is that it sets a fundamental limitation on a deter- 
mination of the Hubble constant from measurements of the 
primary CMB anisotropies alone. This arises because the 
Hubble constant is a secondary parameter, determined from 
the constraint equation (^). The error in the determination 
of h will be dominated by the error in uja, which, as we have 
seen above, is essentially unconstrained because of the ge- 
ometrical degeneracy. The value of h determined from the 
CMB anisotropies is therefore degenerate with the value of 
Qa- For satellite experiments such as MAP and Planck, we 
can assume to a good approximation that the parameters 
ujb, uJc. and ujd are determined exactly. Equation (Q) then 



© 0000 RAS, MNRAS 000, 000-000 



8 G. Efstathiou and J.R. Bond 




I I 

Figure 5. Derivatives of Ci with respect to the seven parameters of the spatially flat target model defined in the text. The derivatives 
with respect to loi, nt and r have been multiplied by the factor indicated in the figure (e.g. the derivative dCi/duJb has been divided by 
a factor of ten). 



provides a relation between h and Qa 
8. 



h = h t /(l- =0a) 



1/2 



(12) 



Here, ht denotes the Hubble constant of the target model. 
Formally, equation (|l^) is valid for small variations around 
the target model parameters, but in practice it is reasonably 
accurate for large variations in Qa (because the degeneracy 
lines in Figure [lib are very nearly straight lines). Equation 
( p^ ) is independent of the details of the experiment. The 
scatter around the relationship (12), however, is set by the 
errors of the density parameters (mainly errors in lud) and 
so depends on the experimental parameters. 

These points are illustrated in Figure ^ which shows the 
results of parameter estimation applied to a set of simulated 
Ci spectra. Here we have generated Ct spectra of our target 
model with errors computed for the OMAP experiment and 
we have determined the cosmological parameters for each 
model by x 2 minimisation, assuming that the theoretical 
spectra of neighbouring points can be approximated by the 
Taylor series expansion 



Ci{si) = C e (s ) + 



( dCt 
V dsi 



Asi. 



(13) 



This procedure is designed by construction to mimic the 
Fisher matrix approach exactly. Figure ^ shows that lu c and 
cjj, are slightly anticorrelated, with errors as given in Table 
The parameters lua and uik are extremely highly corre- 



3.3 



lated along the degeneracy line denned by equation (10b) 



The scatter about this line is determined by the errors in 
ujd- The parameters h and Qa are also strongly correlated, 
following equation Jl2| ) to high accuracy. The final plot in 
Figure W shows the parameters fiA plotted against fife. This 
plot may seem a little strange, because it shows a clump 
of points at Ha ~ 0.8 and fife ~ 0.13, well away from the 



target model parameters fiA = fife = 0. This is because the 
errors in lua and u)d are symmetrical by construction and 
we determine fiA and fife from lua, ud and the Hubble con- 
stant as computed from the constraint relation of equation 
(7). Since the error in lua is large, the points in the Qa~ 
fife plane must, by construction, cluster around the values 
fi A « 7/8, fife « 1/8. 

In summary, the geometrical degeneracy leads to a fun- 
damental indeterminacy in the parameters fiA and fife. For 
spatially fiat target models, this leads to a near indeter- 
minacy of fiA and an error on fife of ~ 0.1 irrespective of 
the accuracy of the CMB experiment. The indeterminacy of 
fiA leads to an indeterminacy of the Hubble constant via 
the constraint relation (7). As a corollary, independent con- 
straints on Ho and fiA from observations other than the 
CMB anisotropics can be used to break the geometrical de- 
generacy. 

The discussion of this Section also shows that tables 
designed to de mon strate the power of various experiments, 
such as Table 3.3, must be interpreted with caution. Evi- 



dently the errors on some cosmological parameters are ex- 
tremely sensitive to the assumed parameter set. For exam- 
ple, restricting to models with zero cosmological constant, 
or curvature, would have a very large effect on the accuracy 
of the determination of Ho- Similar remarks apply to other 
parameters and to near degeneracies other than the geomet- 
rical degeneracy described here. This is the topic of the next 
Section. 



EIGENVECTORS OF THE FISHER MATRIX 
AND PARAMETER DEGENERACIES 
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4.1 Motivation 

The geometrical degeneracy described in the previous Sec- 
tion represents an extreme case of a near exact degener- 
acy between cosmological parameters. In this Section, we 
will investigate partial degeneracies among other parame- 
ters showing, for example, how uncertainties in the param- 
eters defining the primordial spectra couple to the cosmo- 
logical densities in baryons and CDM. The results of the 
previous section show that the geometrical degeneracy is so 
nearly exact that very little is lost in analysing degeneracies 
among other parameters by restricting to a spatially flat 
background universe. Such a restriction merely removes the 
geometrical degeneracy from consideration without affecting 
any other aspect of the analysis. In this Section, we there- 
fore analyse the parameter errors for t he same spatially flat 
target model as defined in Section 3.1 but restricted to a set 
of 7 parameters, n s , rit, Q, r, uit, <^c and u>\. The derivatives 
of Ct with respect to these seven parameters are plotted in 
Figure |E| 

4.2 Principal components 

The Fisher matrix F is a symmetric n x n matrix and so 
can be reduced to diagonal form, 

t 



F 



UAW 



A = diag(Ai, A 2 , A„) . 



(14) 



where U is the matrix in which the m'th row is the eigenvec- 
tor u m corresponding to the eigenvalue A m . The eigenvectors 
are assumed ordered so that Ai > A2 > A3... > A m . From 
our original variables s we can therefore construct a set of 
new variables X that are orthogonal to each other, 

X = U*b, s t Fs = X t AX. (15) 

Thus, for a given experiment and set of parameters s;, we 
can construct a set of variables Xi that are orthogonal lin- 
ear combinations of the original variables; X\ is the most 
accurately determined parameter, X2 is the next most ac- 
curately determined parameter and so on. We call the vari- 
ables X the principal components of the experiment, and 
refer to Xi, Xe etc as 'low-order' components because they 
are poorly constrained (smallest eigenvalues) and X\, X2, 
etc as 'high-order' components. 

We have computed the eigenvalues and principal com- 
ponents for each set of experimental p ara meters listed in 
Table ^. The results are given in Table 



4.1 



which lists the 

components of the projection vector relating the princi- 
pal components t o th e original variables Si. (Note that the 
are computed assuming that the vari- 
have zero mean.) The derivatives of Ce with respect 



numbers in Table 4.1 
ables s 



to the principal components for each experiment are plotted 
in Figure ^| We make the following points concerning this 
analysis: 

[1] As in many aspects of multivariate analysis, certain re- 
sults depend on the scaling of the parameters. The princi- 
pal components in Table JO] are not unique and depend on 
the relative scalings of the variables Sj, as well as on the 
functional form of the variables {e.g. whether we use log U\, 
instead of ujb, or Q 2 instead of Q). One must therefore be 
cautious in assigning a physical interpretation to the prin- 
cipal components. The point of view adopted here (in line 
with most applications of principal component analysis, see 



e.g. Kendall 1975) is that the principal components provide 
a computational tool for assessing whether a set of observa- 
tional points Si he within a sub-space of the n-dimensional 
hyperspace, i = 1, . . . n. (i.e. whether there are degeneracies 
among the physical parameters which are unresolved by the 
experiment) . 

[2] The standard deviation of the principal component Xi 

— 1/2 

is equal to A^ . Thus, the eigenvalues have sometimes 
been used to indicate the power of a CM B experiment (e.g. 
BET97). For example, from Table fll] we see that OMAP 
determines three principal components to an accuracy of 
< 1CP 3 , while CMAP and Planck determine 4 and 5 com- 
ponents respectively to this accuracy or better. While this 
is true, what is important is how many physical (i.e. cosmo- 
logically interesting) parameters can be determined to high 
accuracy rather than the number of linear combinations like 
principal components that can be determined to high accu- 
racy. 

[3] The derivatives of Ce with respect to the physical vari- 
ables plotted in Figure [| can be broadly grouped into two 
categories: those that depend on inflationary variables, n s , 
rit and r that are largest at low multipoles, and those that 
depend on physical densities, u)b, w c , to a, that have signifi- 
cant amplitudes to I ~ 2000 — 3000. Cosmic variance sets a 
fundamental limit to the accuracy with which the variables 
n s , rit and r can determined from observations of the CMB, 
whereas the cosmic densities can be determined to high pre- 
cision if Ce is measured accurately at high multipoles. Thus, 
the power of an experiment can be gauged by the extent 
to which the principal components (particularly those with 
the largest variances) mix inflationary variables with cos- 
mic densities. Qualitatively, this can be assessed from Fig- 
ure |. For OMAP and CMAP, the C e derivatives with re- 
spect to the low-order principal components are oscillatory 
and have a significant amplitude at multipoles £ ~ 100 - 
2000. This tells us immediately that the low-order principal 
components involve the cosmic densities at some level and 
hence that they will be partially degenerate with inflation- 
ary variables, which is obviously unsatisfactory if we want 
to measure the cosmic densities to high accuracy. 
[4] The principal components for a Planck-type exp erim ent 
are nearly equivalent to physical variables. In Table 4.1, we 



have placed an asterisk next to each entry with an abso- 
lute value greater than 0.5, so that the reader can see easily 
which physical variables contribute strongly to each princi- 
pal component. For Planck, Xi, X%, A3, X4 and A5 couple 
strongly to Ub, uj c , Q, n s and ua respectively. The two low- 
est order principal components are made up predominantly 
of the two tensor parameters r and n± . Nevertheless, the de- 
generacies that remain dominate the errors of the cosmolog- 
ical densities. These remaining degeneracies are, in a sense, 
fundamental, because they cannot be removed by making 
more accurate observations of the linear CMB temperature 
anisotropies. The results for Planck shown in Table JO] and 
Figure n are very close to those for an experiment limited 
by cosmic variance alone. Accurate CMB polarization mea- 
surements on large scales could, however, constrain the ten- 
sor component and so improve the accuracies of the cosmic 
densities. Section demonstrates how the accuracies of the 
cosmic densities uib and ui c are affected by constraints on the 
tensor component. 

In the rest of this Section, we investigate the degen- 
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Table 3: Principal components 



OMAP errors 







2.27E-4 


0.998* 


3.14E-3 


5.363E-2 


4.46E-3 


2.727E-2 


1.91E-2 


2.722E-2 


4.32E-2 


-6.283E-3 


3.18E-1 


3.315E-3 


1.06 


6.271E-4 



U) c 




-1.417E-2 


-1.103E-2 


-0.284 


3.051E-3 


0.886* 


0.321 


0.135 


0.214 


-0.312 


0.821* 


-0.134 


0.416 


-1.768E-2 


5.901E-2 



Q 


n s 


-6.376E-2 


4.919E-3 


0.871* 


-0.397 


1.513E-2 


-0.294 


0.424 


0.829* 


-0.183 


-0.191 


1.798E-2 


0.177 


4.918E-3 


2.753E-2 



r 


nt 


-2.972E-4 


-2.599E-4 


1.258E-2 


1.188E-3 


2.910E-2 


2.633E-2 


-0.198 


-0.167 


0.300 


-0.261 


0.565* 


0.676* 


0.742* 


-0.667* 



CMAP errors 





u b 


1.89E-4 


0.998* 


2.98E-3 


4.883E-2 


3.82E-3 


3.571E-2 


1.22E-2 


2.471E-2 


3.39E-2 


2.270E-3 


2.46E-1 


3.025E-3 


0.997 


-2.923E-4 







-1.592E-2 


-1.113E-2 


-0.366 


-3.678E-2 


0.870* 


0.314 


0.149 


-6.142E-2 


-0.261 


0.849* 


-0.136 


0.417 


1.032E-2 


-3.328E-2 



Q 


n s 


-6.241E-2 


3.257E-3 


0.862* 


-0.344 


0.279 


-0.252 


0.415 


0.886* 


-4.378E-2 


5.800E-2 


1.385E-2 


0.171 


-1.769E-3 


-1.260E-2 



r 


nt 


-2.052E-4 


-1.795E-4 


1.005E-2 


9.545E-3 


2.238E-2 


2.0370E-2 


-9.596E-2 


-8.279E-2 


-0.345 


-0.295 


0.616* 


0.631* 


-0.701* 


0.712* 



PLANCK errors 







4.09E-5 


0.999* 


5.29E-4 


3.195E-2 


7.42E-4 


-3.704E-2 


4.32E-3 


-1.172E-2 


5.25E-2 


7.997E-3 


5.54E-2 


-5.331E-4 


0.75 


-1.424E-5 



U) c 




-3.110E-2 


-1.707E-2 


0.958* 


0.274 


-5.976E-2 


-4.110E-2 


-3.611E-2 


7.679E-2 


-0.275 


0.957* 


-4.795E-4 


-2.665E-3 


1.003E-4 


-4.164E-4 



Q 


n s 


3.253E-2 


1.847E-2 


6.588E-2 


2.740E-2 


0.977* 


0.196 


-0.195 


0.976* 


3.922E-2 


-7.733E-2 


-1.017E-2 


4.461E-2 


-2.579E-4 


9.970E-4 



r 


nt 


-1.088E-5 


-9.348E-6 


2.852E-4 


2.541E-4 


8.106E-4 


7.508E-4 


-3.396-2 


-3.025E-2 


3.716E-3 


3.843E-3 


0.760* 


0.648* 


-0.649* 


0.761* 



Table 1. Components of the projection vector relating the original variables §i (with means subtracted, s; = Si — (si)) to the principal 
components X. The eigenvalues are listed in the first column. We have marked each component of U* that is greater than 0.5 with an 
asterisk to show which physical variables contribute significantly to each principal component. 



eracies among cosmological parameters as denned by the 
low-order principal com ponents. We conclude this Section 
by listing in Table 4.2, under the column headings 'all pc', 
the la errors in the cosmological parameters derived from 
the Fisher matrix for the spatially flat case. 

For most of the parameters, the error estimates are cl ose 
to those listed in Table B.3. The errors on lua in Table 4.2 



4.3 Degeneracies among cosmic parameters 



are obviously much smaller than those given in Table i.3 be- 
cause we are restricting to a spatially flat background model. 
The erro rs on uja will be discussed in greater detail in Sec- 
tion 4.4. The largest differences between Tables 3.3 and 4.2 



are for the errors in the tensor parameters r and nt- The 
errors on these parameters are very large, so the first or- 
der expansion of CV, on which the Fisher matrix expression 
(4) is based, cannot necessarily be expected to give accurate 
answers. 



The columns in Table 



4.2 



labelled 'lpc', '2pc' and '3pc', list 
the variances of the cosmological parameters that arise if we 
include in the Fisher matrix analysis only the lowest order 
principal component (Xr, column labelled 'lpc'), the two 
and three lowest order principal components (Xg and Xr in 
the column labelled '2pc', and X$, X& and Xr in the column 
labelled '3pc'). For CMAP and OMAP, the two lowest order 
principal components account for most of the variance of all 
the cosmological parameters listed in the table. Since we have 
specialized to a spatially flat universe, the geometrical de- 
generacy described in Section [| has been explicitly removed 
from the analysis. Neve rthe less, to a good approximation, 
the parameters of Table [4.2| lie on a two dimensional surface 
within the 7-dimensional parameter space. It is the orienta- 
tion of this surface with respect to the axes defined by the 
parameters that accounts for the most of their variance. 

This is illustrated in Figure [?]. In Figure 0a, we have gen- 
erated points in the 7-dimensional space of our cosmological 
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Figure 6. Derivatives of with respect to the principal components listed in Table 4.1 for the OMAP, CMAP and Planck experiments. 
The principal components have been ordered so that X\ is the most accurately determined (i.e. has the largest eigenvalue), X2 is the 
next most accurately determined and so on. The derivatives of Ci with respect to X\ and X5-X7 have been multiplied by the factors 
indicated in each figure. 



para mete rs by using all principal components as given in 
Table 4.1 for the OMAP experiment. The principal compo- 



nents are assumed to be Gaussian distributed with disper- 
sion equal to l/VAi- We have made one significant change 
to the Fisher matrix analysis of Section Q in that we have 
applied a constraint that the densities uji and amplitude r 
must be positive. This is evidently physically reasonable for 
the uji, required for r, and aids in the comparison with the 
exact calculations shown in Figure^, to be described below. 



We show correlations between various pairs of parameters, 
but note in particular the very strong correlations between 
wa and ui c and between h and S1a- The points in these plots 
lie tightly on lines, similar to the situation shown in Figure 
^, which illustrated the geometrical degeneracy. As we will 
show below, the strong correlations between uja and uj c and 
h and Qk have a similar physical explanation to the geomet- 
rical degeneracy, though the degeneracy is not exact in this 
case. 
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(a) all principal components Fisher matrix 



(b) two principal components Fisher matrix 
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(c) two principal components exact 
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Figure 7. Monte-Carlo simulations of parameter estimation as described in the text showing correlations between various pairs of 
parameters. We assume the experimental parameters for OMAP as given in Table [| In panel (a) we generate simulations using all 
principal components, reproducing the variances and correlations compute d fro m the Fisher matrix. In panel (b) we generate simulations 
using only the two lowest order principal components Xq and X-j (Table |3.3[ ). In panel (c) we show the results of an exact likelihood 
analysis applied to a two-dimensional grid of models in the variables X r an d Xi f see figure |j and the discussion in the text). The solid 
lines in the u;a — w c and h — Qa panels are computed from equations ( ^3a| ) and ( |231: ) . The das hed line in the uic—uif, plane shows the 
degeneracy line given by a constant height for the first Doppler peak as described in Section |4,5} 
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Table 4: lcr errors in estimates of cosmological 
parameters (spatially flat universe) 



OMAP 



cannot allow negative values of u>b, u> c and r and (in this 
example) we do not allow negative values of wa- The cross 
in Figure Ha shows the ±la errors on Xe and X7 computed 



Porn Tt~t nf or 


all pc 


lpc 


2pc 


3pc 


5u> b /u b 


0.11 


0.056 


0.098 


0.10 


5u) c /u) c 


0.21 


0.080 


0.19 


0.20 


SbJA 


0.15 


0.064 


0.14 


0.15 


SQ 


0.014 


0.0052 


0.0079 


0.011 


Sr 


0.81 


0.81 


0.81 


0.81 


Sn s 


0.066 


0.030 


0.062 


0.063 


Sn t 


0.745 


0.73 


0.74 


0.74 






CMAP 




l.V y. jji yy. 

1 al CLI11(_- LL 1 


all 

a pc 


lpc 


2pc 


pc 


5uj b /uj b 


0.072 


0.024 


0.063 


0.064 




0.15 


0.044 


0.15 


0.15 


Sua 


0.11 


0.034 


0.11 


0.11 


SQ 


0.0071 


0.0018 


0.0038 


0.0041 


Sr 


0.70 


0.70 


0.70 


0.70 


Stl s 




0.0013 


0.044 


044 


Snt 


0.71 


0.71 


0.71 


0.71 






PLANCK 




Parameter 


all pc 


lpc 


2pc 


3pc 


5u) b /u> b 


0.0068 


0.0008 


0.0024 


0.0040 


5u) c /u) c 


0.0063 


0.0003 


0.0003 


0.0060 


8w A 


0.0049 


0.0003 


0.0003 


0.0049 


SQ 


0.0013 


0.0002 


0.0006 


0.0006 


Sr 


0.49 


0.49 


0.49 


0.49 


Sn s 


0.0050 


0.0008 


0.0026 


0.0026 


Snt 


0.57 


0.57 


0.57 


0.57 



Figure 0b is constructed in exactly the same way as 
Figure up,, except that we include only the two lowest order 
principal components. As expected from the entries in Ta- 
ble 4.2, by incorporating only the two lowest order principal 



components, we reproduce the degeneracy lines in Figure |7p 
which account for most of the variance of the physical pa- 
rameters. The higher order principal components contribute 
to the scatter around these degeneracy lines but do not con- 
tribute much to the total variance. 



Table 4.2 shows that the errors on some physical vari- 
ables are large (e.g. uia, r, nt). As a consequence, the vari- 
ances of the low order principal components are also large 
and so this raises a question concerning the validity of the 
Fisher matrix approach to parameter errors, since this relies 
on small variations Asi of the parameters around those of 
the target model. To test the validity of the Fisher matrix 
approach, we have computed an exact version of Figure 0b 
using maximum likelihood. We define a 400 x 400 grid in the 
principal components Xe and X7 and compute scalar and 
tensor power spectra using the CMBFAST code. The ranges 
of Xe and X7 are illustrated in Figure ^a. The contours in 
Figure []a show 1, 2 and 3cr contours for the likelihood ratio 
computed from equation ([)]) with £ ma x = 1000. The jagged 
lower boundary of the contours in the figures is imposed 
by positivity constraints on the physical parameters, i.e. we 



from the eigenvalues listed in Table 4.1; note that the scales 
of Xe and X7 in this figure have been chosen so that the er- 
ror contours computed from the Fisher matrix are circular. 
Qualitatively, the likelihood errors are of about the same 
size as those computed from the eigenvalues, but the dis- 
tributions computed from the likelihood function are more 
sharply peaked than a bivariate Gaussian distribution. This, 
together with the positivity constraints which further nar- 
rows the distribution, means that the errors in Xe and X7 
computed from the likelihood ratios are smaller than those 
computed from the eigenvalues of the Fisher matrix. 

Figure ^b shows a series of Monte-Carlo simulations. 
For each point plotted in the figure we generated a sim- 
ulated power spectrum Ci assuming Gaussian errors with 
variance given by equation ( pa[ ) with the parameters of the 
OMAP experiment listed in Table 3.3. For each simulated 
Ce spectrum, we computed the maximum likelihood values 
of Xq and X7 by comparing with the grid of numerical mod- 
els computed with CMBFAST. (In fact, we reduce the effect 
of the finite grid using bi-linear interpolation to compute 
the maximum likelihood values of Xe and X7.) Each point 
in Figure ||d is plotted at the best fitting value of Xe and 
X7, and the distribution of points follows the likelihood con- 
tours extremely well showing that our computations are self 
consistent. From the values of Xe and X7, we compute the 
physical parameters Si and so we can plot an analogue of Fig- 
ure ^b using parameters computed from an exact likelihood 
approach. This is shown in Figure 0o. The general trends 
between the variables are similar to those shown in Figures 
(7^ and 7b, but the ranges are significantly narrower showing 
that the Fisher matrix approach overestimates the errors of 
the physical parameters. The Fisher matrix approach does 
give qualitatively correct results (which is impressive given 
the large errors in some of the physical variables) but evi- 
dently does not give precise error estimates. 

Furthermore, since much of the variance in all of the 
physical parameters comes from variations of the two low- 
est order principal components, the Fisher matrix approach 
will overestimate the errors in apparently well determined 
parameters such as u) b , since these are dominated by the 
variances of the low-order principal components. This is also 
illustrated by the entries in Table 4.2 which is a tabular 
version of Figure 0. The second and third columns list the 
variances of the physical parameters (including h and SI a) 
computed from all principal components and the two low- 
est order principal components; these numbers d iffer from 
the corresponding numbers for OMAP in Table 4.2, only 
because we have applied positivity constraints on r and the 
densities uja, lj c and uj b . The fourth column shows the lcr 
errors in the physical parameters computed from the Monte- 
Carlo simulations plotted in Figure [7)3. Table 4.5 shows that 
the Fisher matrix overestimates the errors in many of these 
parameters by a factor of two or so. 

The example shown in Figure [?]c and Table 4.3 shows 
how the principal component analysis can be used to sim- 
plify a likelihood calculation of the errors in cosmological pa- 
rameters, so improving on the Fisher matrix approximation. 
In general, the cosmological parameters defining a model 
are highly degenerate with respect to each other and so the 
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Figure 8. Likelihood ratio contours in the X§-Xf plane for the 
OMAP experiment (1, 2 and 3cr contours computed as described 
in the caption to Figure 3). The cross in Figure Ha shows the 
la errors derived from th e Fi sher matrix (computed from the 
eigenvalues listed in Table |4.l| ). The four filled circles show the 
locations of the four models plotted in Figures [ll] and [12] below. 
Figure |^b shows the maximum likelihood estimates of Xq and X7 
for the Monte-Carlo models described in the text. These values 
of X§ and X-j were used to compute the values of the physical 
parameters plotted in Figure f^c. 

variances of the s; are dominated by a small number of p rin- 
cipal components. In the examples given in Table 4.2, the 
two lowest order principal components account for most of 
the variance of the physical parameters. One can therefore 
use the principal components of the Fisher matrix to reduce 
the dimensionality of the parameter set to say two or three 
dimensions. It is then feasible to evaluate the CVs numeri- 
cally on a grid in this reduced parameter space using a fast 
code such as CMBFAST. The numerically computed power 
spectra can then be used to calculate the errors on the phys- 
ical parameters via a likelihood analysis in the reduced di- 
mensional space. This procedure, although time consuming, 
is much faster than attempting to calculate the CVs numer- 
ically on a grid in all of the physical parameters denning the 
model (7 dimensions in the simplified example given above, 



Table 5: 1<t errors in estimates of cosmological 
parameters: comparison of Fisher matrix with 
exact likelihood including positivity constraints. 



OMAP 



Fisher 



Likelihood 



Parameter 


all 


2pc 


2pc 




0.082 


0.061 


0.032 


SuJc /Wc 


0.13 


0.12 


0.071 


SOJA 


0.094 


0.091 


0.052 


SQ 


0.013 


0.0047 


0.0022 


Sr 


0.55 


0.55 


0.24 


Sn s 


0.045 


0.040 


0.022 


Sn t 


0.58 


0.58 


0.32 


8h/h 


0.11 


0.11 


0.064 


sn A 


0.19 


0.18 


0.15 



10 or more dimensions in the more realistic models described 
by BET97 and ZSS97). 

Another consequence of the strong degeneracies among 
physical parameters is that the variance estimates of some 
parameters are extremely sensitive to imposed constraints. 
For example, if we restrict the range of possible models to 
those that have negligible tensor mode contribution (i.e. we 
constrain the parameter r to be zero) then this has a large 
effect on the lowest order principal components and hence 
on the variances of most of the physical param eter s. This 
point will be described in more detail in Section 4.6. 



4.4 



Location of Doppler peaks and degeneracies 
for spatially flat models 



In Section 3.2 we analysed the geometrical degeneracy be- 
tween the parameters LOk and uja- In this section, we analyse 
a similar degeneracy for spatially flat models, i.e. we restrict 
to u>k = and derive the relationship between uja and u> c (or 
ujm) that leads to a nearly identical location of the Doppler 
peaks as a function of multipole. However, since the Doppler 
peak structure (heights and relative amplitudes of the peaks) 
depends on the matter content of the Universe, the resulting 
degeneracy is not exact. As we will see below, the locations 
of the Doppler peaks can explain the strong correlations in 
the lo a -uj c and h-flA plots shown in Figure [?]a. However, a 
precise experiment probing high multipoles such as Pl anc k 
can break these degeneracies to high accuracy (c/ Table 4.2). 



The comoving wavenumber of the m'th Doppler maxi- 
mum or minimum is given approximately by 

k m r a (a r ) — m-K (16) 

(Hu and Sugiyama 1995), where r„ is t he sound horizon at 
recombination as defined in Section 
is given by, 

da 



3.2. The sound horizon 



V3 H nT 

1/2 



V3 HoSl n 



/2 



(a + a equ y/i (l + i?) 1 / 2 

I s (LU m ,LU b ) . 



(17) 



The expression on the right hand side of this equation defines 
the integral I s (oj m , ojt) which is used in equation (E3) below. 
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The parameter a equ is the scale factor at the time that mat- 
ter and radiation (including massless neutrinos) have equal 
density, and R = (3p(,/4p 7 ). Thus, 



-l 



24185 



1.6813 

1 + 7J W 



R = 30496u b a . (18) 



where r\ v denotes the relative densities of massless neutrinos 
and photons, p v = r) v p 7 , and is equal to 0.6813 for three 
massless neutrino flavours. Evaluating the integral in equa- 
tion ([[?]) gives, 

4000 alii 



,1/2 (1+^)V2 



In 



(1 + Riz,.)) 1 ' 2 + (R(z r ) + R equ ) 1/2 



1 + vArT 



Mpc. (19) 



We define the redshift of recombination z r as the redshift at 
which the optical depth to Thomson scattering is unity, as- 
suming the standard recombination history (Peebles 1968). 
A useful fitting formula for z r is given by Hu and Sugiyama 
(1996), 



Zr 

91 

f)2 



1048[1 + 0.00124^ 
0.0783a;, 



[1 + giu% 



(20) 



a238 [l + 39.5a; 6 a763 ]- 1 



0.560[1 +21.1w 6 181 ] 



The locations of the Doppler peaks are therefore given by 
■ dA ^l, (21a) 



l m « mir- 



where cLa is the angular diameter distance to the last scat- 
tering surface 



cLa 



sinh(w£ /2 2/) , 



(21b) 



and y is given by equation (pb[). In fact, equation (21a) 
is approximate because the exact numerical coefficient re- 
lating £ m and o\a /r s depends on the projection of the 
three-dimensional temperature power spectrum to a two- 
dimensional angular power spectrum. The exact relation de- 
pends therefore on the shape of the primordial power spec- 
trum and on the Doppler peak number m (see e.g. Hu and 
White 1996). From the numerical computations of C'e we 
find that the location of the maximum of the first Doppler 
peak for a scale-invariant primordial spectrum of scalar fluc- 
tuations is well approximated by the expression 



0.746ttV3(1 



U/2 



V* sinh[^ 2 y ] 



(22) 



Equation (|2^) is valid for negatively curved models, 
though we specialize to spatially flat models in this section. 
To derive degeneracy lines for spatially flat models, we take 
the limit w& — *• in equation (^). For a fixed value of cot, 
the condition of fixed in equation ( p^ ) then leads to a 
relation between loa and u c . Figure ^| shows the degeneracy 
lines for a spatially flat universe with the baryon density of 
our target model (ub = 0.0125). 

As in Section 3.3, for small variations of the parameters 
around those of the target model, we can derive the de- 
generacy lines by differentiating Id- We do not expect the 
degeneracy lines derived in this way to be exact because the 
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Figure 9. The position of the first Doppler peak, ip, for a spa- 
tially fiat universe {uik = 0) and assuming tjj, = 0.0125. The line s 
show contours of constant £jj , computed from equation ( 23a 231: ) , 
plotted as a function of ua and lo c . The value of Id is given next 
to each contour. 



degeneracies involve the entire Doppler peak structure ac- 
cessible to the experiment, not just the position of the first 
Doppler peak. However, provided the baryon density is de- 
termined to high precision from the experiment (i.e. ut, in 
equation [22] is held fixed), we would expect that the posi- 
tions of Doppler peaks would define the degeneracy lines in 
the io c -lja plane quite accurately. These are given by 

, _ {d£ D /dcu A ) t 



LUc = (Wo)( + bulA, 



(23a) 



where the subscript t denotes that quantities evaluated as- 
suming the parameters of the target model. For our spa- 
tially flat target model, b = —0.368. As in the derivation of 
equation (^), the Hubble constant is degenerate with Q\ 
according to 

ft= (i-d + W a - (23b) 

The solid lines in the lua-u> c a nd h -Q.A panels in Figure [?] 



are computed from equations (23a 23b). As expected, they 
give a good approximation to the degeneracy lines, but the 
agreement is not exact. 

4.5 Height of the first Doppler peak and 
parameter degeneracies 

In analogy with the location of the Doppler peaks, the height 
of the first Doppler peak can provide some further insight 
into parameter degeneracies. The scalar component of the 
power spectrum can be roughly approximated by the fol- 
lowing expression 



e(£+i)Ct e{e + i)r(e+ 



r( 



6C*2 



T{£- 



(5-n s ) ■ 
2 



r( 



(3+n s ) n 
2 
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Figure 10. Simulations of parameter estimation for an experi- 
ment sampling only the first Doppler peak. We assume that the 
experiment is cosmic variance limited over the multipole range 
2 < £ < 300. The solid lines in the u)^-u) c and h—Q\ pl anes show 
lines of constant peak location as given by equations (23a,) 23b 
The dashed lines in the ui c — uii, planes and in the n s — ujf, plane show 
lines of constant Doppler peak height as given by equation 24a . 



A(uj bl uj c ,n a ) 
+ , .,, exp 



{i-i D f 



2(Al D ) 



where 



Air 



0A21 



D. 



(24a) 



(24b) 



]nA(u!b,U!c,n s ) ~ 4.5(n s — 1) +ai + <i2UJ c + astuc (j24c) 



ai = 2.376, a 2 = 3.681, a 3 = -5.408, 
a 4 = -54.262, a 5 = 18.909, a 6 = 15.384. 



(24d) 



and lp is given by equation ( |22[ ) The first term in equation 
(24a) is the usual Sachs-Wolfe expression for a power-law 
fluctuation spectrum (see e.g. Bond 1996). The second term 
is a fit to the first Doppler peak determined from a grid of 
CMBFAST computations of C t . 

The main parameter degeneracies for an experiment 
sampling only the multipole range £ ,$ 300-400 will be de- 
termined by the properties of the first Doppler peak. Thus, 
the Doppler peak position £d will determine the degener- 
acy between u>a, t^c and w& (as described in the previous 
section) and the Doppler peak height will determine the de- 
generacies among ujb, u c and n 3 . This is illustrated in Figure 
|ic| , which shows a similar plot to Figure 0a, but for an ex- 
periment that is cosmic variance limited for the multipole 
range 2 < £ < 300. The conditions of constant peak height 
and location provide a good description of the parameter 
degeneracies in this case. 

The situation is more complicated for a MAP or Planck 
type experiment because they are sensitive to the subsidiary 
Doppler peak structure. This is illustrated by the uj c -u)b 
scatter plot in Figure ^a. To a reasonable approximation, 
the subsidiary Doppler peaks pin down the value of u>b to 
high accuracy and so the first principal component is almost 



vertical in the uj c -ui b plane (in fact it is tilted to a slightly 
negative slope to reduce the scatter in the Doppler peak lo- 
cations). In this case, the condition of constant first Doppler 
peak height (plotted as the dashed line) determines the scat- 
ter in the orthogonal direction to that denned by the first 
principal component. 



4.6 Constraints on the Tensor Component 

Figure [H] shows the temperature power spectra for the four 



models indicated by the dots in Figure |8ja. These models 
lie on the 2er contour in the plane defined by the lowest or- 
der principal components. An experiment with the OMAP 
parameters can therefore only marginally distinguish these 
models from our target model. Figure til is designed to il- 
lustrate the following points: 

[1] The parameters for these models are given in the up- 
per panels. Notice that the the tensor parameters nt and 
r vary by large amounts and that the most obvious differ- 
ences with the power spectrum of the target model are at 
low multipoles. This is expected because the poorly deter- 
mined tensor parameters carry hi gh w eight in the lowest two 
principal components (c.f. Table 4.2). 



[2] The dotted lines in the upper and lower panels show the 
power spectra computed from the first derivatives according 
to equation (^). These differ significantly from the numer- 
ical computations (shown as the dashed lines) because the 
variations in nt and r, and of some other parameters, are 
so large that the first order Taylor expansion breaks down. 
This is why the exact likelihood analysis shown in Figure ^c 
differs from the Fisher matrix approximation used to derive 
Figure 0b. The exact analysis leads to tighter constraints 
on the physical parameters because the first order approxi- 
mation tends to underestimate the differences in the power 
spectra. 

[3] Throughout this paper, we have allowed the parameters 
n s , nt and r to vary independently of each other. In other 
words, we have not imposed any additional constraints on 
the physical mechanism which gave rise to the fluctuations, 
other than that they are Gaussian characterised by power- 
law spectra. However, wide classes of inflationary models 
lead to constraints between the scalar and tensor compo- 
nent (see the comprehensive review by Lidsey etal 1997 and 
references therein). Such constraints can severely limit the 
allowed area in the plane denned by the lowest two principal 
components, and since these account for such a large pro- 
portion of the variance of all the physical parameters, the 
parameter errors will be extremely sensitive to these con- 
straints. 

We will investigate this last point in the remainder of 
this Section, using as examples two limiting constraints. It is 
useful first to review some key results concerning the scalar 
and tensor perturbations predicted by inflationary models 
to motivate these approximations. 

Tensor perturbations, arising from quantum noise in the 
gravitational wave degrees of freedom, obey a perturbation 
equation like that of a free massless field, whereas the equa- 
tion for scalar curvature fluctuations have extra terms deter- 
mined by the shape of the scalar field potential V((j>), where 
<j> is the inflaton field that drives inflation. The relationship 
between the two is therefore model dependent, and can be 
quite complex if V(<f>) is, e.g., if inflation involves more than 
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Figure 11. Temperature power spectra for the four models plotted as the filled circles in Figure ga. For these models the two least 
well determined principle components Xq and Xj lie on the 2<r contour in Figure ||a. The physical parameters of these models are given 
in each panel. The solid line in the upper panel in each figure shows the power spectrum of the target model; the dashed line shows 
the exact computation of Ct using CMBFAST for the parameters defined by X% and X7. The dashed line shows Ct inferred from the 
derivatives of Ct (equation ha). The lower panels show the residuals of these curves with respect to the target model; dashed lines show 
the residuals for the numerical computations of Ct and dotted lines show the residuals of the derivative approximation. The lines labelled 
OM, CM and P show the la errors computed from equation (pa|) for the parameters of OMAP, CMAP and Planck given in Table Li 



one dynamically important scalar degree of freedom. The es- 
sential features of inflation for one degree of freedom can be 
understood using a Hamilton- Jacobi formulation (Salopek 
and Bond 1990) in which the Hubble parameter H(<f>) is 
treated as a function related to V(<f>) by 



H 1 



3rrip 



4tt 



m ^> +V(<t>) 



(25) 



where m p is the Planck mass. The amplitudes of the post- 
inflation scalar and tensor fluctuation spectra, Ps and Pt, 
can be expressed in terms of the value of H when a comoving 
wavenumber k crosses the Hubble radius k — Ha at time rt : 



Ps(k) 



1 



1 + q 



H 2 (r k ) 2us 



(26a) 



P T(k) = ^!il*) e 2« 4 



(26b) 



The deceleration parameter q is related to the slope of H by 

2 



(1 + 9) = 



ml ( d\nH \ ' 



47r V 



(27) 



(e.g., Bond 1996). The parameters u s and ut (e.g., Bond 
1994) turn out to be small for power-law like potentials 
(but testing whether they are non-zero is an important goal 
of CMB measurements). Thus the ratio of tensor to scalar 
power is largely determined by (1 + q), which measures the 
departure from pure exponential expansion during inflation. 

Taking logarithmic derivatives of the spectra Ps and Pt 
with respect to wavenumber k gives the scalar and tensor 
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spectral indices n 3 and nt introduced in Section 

,2 «2, 



n s = 1 + 2 1 



1 m; a 2 injj- 

A2~ +£s ' 



g y g 27r <9<j!> 2 



n t = 2 ( 1+ - ) +e t , 



(28a) 
(28b) 



where the terms e s and et are small corrections arising from 
derivatives of u s and ut. Thus the difference n s — f — nt is de- 
termined by the change in q over the k range appropriate for 
large scale structure and the CMB, which turns out to cor- 
respond to quite a narrow section of the inflaton potential. 
The familiar 'slow rollover' approximation assumes q ~ — 1 
in equation (pBj), giving H 2 « 8 ttV / (3m 2 ,) , but not in equa- 
tion (26a) nor in equations ( 28a , |28b ) . Indeed for inflation to 
end, the universe must enter a deceleration phase, so realis- 
tic potentials must give departures from pure scale-invariant 
spectra, possibly large but more plausibly small. 

The ratio of the scalar and tensor components can also 
be computed in terms of the inflation parameters defined 
above; e.g., if n„ and nt are assumed to be constant, 



r = 13.7(1 + q)e 



2(u t — «s)g— 0.15n t g — l.l(m s 



-6.8n t 



e -0.15nt e -l.l(n 3 -l-n t ) 



'"^29) 
(30) 



(Bond 1996). r also explicitly depends upon parameters such 
as and uj^- 

There are two special cases we use to illustrate con- 
straints on the tensor component. The first is nearly uniform 
acceleration, (so d 2 In H/d(f> 2 ~ 0), 

e.s ^-^lm a) 



-6.8 



n t 



(l-n t /2) 



(3 - n 3 )/2 v 



Pure power law expansions have q precisely constant 
whereas chaotic inflation models with smoothly varying 
power law potentials have q nearly uniform, but the devia- 
tions from n s — 1 ~ n t are important to include if n t is small. 
Equation (31a), often approximated by r ~ 7(1 — n 3 ), has 



been discussed by many authors (see e.g. Lidsey et al. 1997, 
and the footnote in Section 3.1). The second constraint we 
consider has 



nt ~ 0, 



0, 



(31b) 



which applies to models such as 'natural inflation' (Adams 
et al. 1993) in which the inflaton begins near the maximum 
of an extremely flat potential, so 1 + q and the tensor mode 
are exponentially suppressed. Since d(l + q)/d(f> cannot then 
be neglected, n„ can be less than unity, possibly significantly 
so. Similar behaviour can arise in inflationary models based 
on supergravity (e.g. Ross and Sarkar 1996). For specific 
classes of models, one can refine the estimates of equations 
(Sis ,31b), for example by relating the corrections u t — u s to 
the spectral i ndic es. 
In Table 



4.6 



we show how the errors in the cosmological 
paramete rs change if we apply the constraints of equation 
(Sis, 31b). The table lists the la errors computed from the 



Fisher matrix including positivity constraints. The columns 
labelled 'no const.' give the results when no constraints are 
applied to the tensor spectrum, i.e. all parameters are al- 
lowed to vary independently. The entries under this heading 
for OMAP are thus identical to the entries in the second col- 



umn of Table 4.3. The columns labe lled V = —7nt show the 
results of applying the constraint (Bla) and those labelled 



'r = 0' show the result of applying the constraint (311) 



As expected (see point 3 above), imposing the con- 
straints on a OMAP-type experiment results in a signif- 
icant reduction of the errors on most parameters. Even 
for a CMAP-type experiment, the constraints lead to large 
changes in the errors of some parameters, e.g. the fractional 
error in uid decreases by about a factor of 3, the fractional 
error in h decreases by a factor of 4 and the error in £Ia de- 
creases by a factor of between 2 and 3. For a Planck-type ex- 
periment, however, imposing the constraints le ads to hardly 
any change in the parameter errors. As Table 4.1 and Fig- 
ure ^| show, for a Planck-type experiment, there is almost a 
decoupling between the highly uncertain tensor parameters 
and the high order principal components. 

The results of Table [L(j show that the errors on cosmo- 
logical parameters can be reduced by making certain specific 
assumptions concerning the inflationary model. For OMAP 
and CMAP-type experiments, the reduction in the errors 
is very significant, but at the expense of overly restricting 
the inflation model. Allowing full freedom in how V(4>) and 
thus H(4>) i s structured can result in strong vari ati ons in q, 
with the consequence that the power ratio 26b to 26a and 
the indices n s (k) and n t (k) could be complex functions of 
the wavenumber. As might be expected, if there is complete 
functional freedom, the errors on the cosmological parame- 
ters become much larger (e.g., Souradeep et al. 1998). Such 
inflation models with many features over the small regime 
of V((j)) that the CMB can probe seem rather baroque, and 
constraining the potential by invoking smoothness "prior 
probabilities" seems perfectly reasonable. For example, a 
simple way that this has been included in the past is to 
allow the spectral indices to "run" with a logarithmic cor- 
rection: n s (k) « n a (k s ) + a \n(k/k s ) (Kosowsky and Turner 
1995, Lidsey et al. 1997). Of course, in presenting the re- 
sults of CMB parameter estimation, the constraints (if any) 
applied to the tensor mode must be spelled out precisely. 

This analysis illustrates also that constraints on the 
tensor component derived from CMB polarization measure- 
ments, could lead to a significant improvement in the esti- 
mates of other parameters, including the baryon and CDM 
densities. 



5 STRENGTH, SIGNIFICANCE AND 
CORRELATED ERRORS IN C t 

The next generation of CMB satellites should lead to esti- 
mates of Ct with errors dominated by cosmic variance, 



AC e 
Ct 



(2£+l)f 3ky 



for multipoles £ ~ 1000 - 2500. At £ ~ 1000, cosmic variance 
limits the la accuracy of a measurement of Ct to about 3 
percent. At first sight, it might seem as though there is little 
point in measuring Ct (and subtracting contaminating fore- 
grounds) to a much higher accuracy than the limit imposed 
by cosmic variance. However, such thinking would be incor- 
rect. The reason is related to the concepts of strength and 
significance in statistics (see Press et al. 1996) and we illus- 
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Table 6: la errors in estimates of cosmological parameters: 
effect of constraining the tensor component. 







OMAP 






CMAP 






PLANCK 




Param. 


no const. 


r = — 7n t 


r = 


no const. 


r = —7th 


r = 


no const. 


r = ~7n t 


r = 


5u b /u) b 


0.082 


0.036 


0.046 


0.052 


0.028 


0.030 


0.0064 


0.0056 


0.0056 


Sw c /w c 


0.13 


0.046 


0.044 


0.097 


0.028 


0.031 


0.0042 


0.0042 


0.0039 


Suj a 


0.094 


0.025 


0.026 


0.069 


0.014 


0.020 


0.0030 


0.0030 


0.0027 


SQ 


0.013 


0.007 


0.009 


0.0066 


0.0047 


0.0050 


0.0013 


0.0010 


0.0011 


5r 


0.55 


0.054 




0.49 


0.043 




0.33 


0.023 




Sn s 


0.045 


0.0077 


0.015 


0.030 


0.0061 


0.0098 


0.0049 


0.0032 


0.0042 


Sn t 


0.58 


0.0077 




0.56 


0.0061 




0.40 


0.0032 




Sh/h 


0.11 


0.032 


0.033 


0.082 


0.020 


0.028 


0.0045 


0.0045 


0.0041 


sn A 


0.19 


0.079 


0.082 


0.16 


0.049 


0.068 


0.012 


0.012 


0.011 



trate the point with the following simple example. Imagine 
that we have N independent random variates, each drawn 
from a Gaussian with a mean and dispersion of unity. From 
these N random variates we can estimate the mean 

with variance 

Var( M ) = 1. 

Thus, even though the variance on each individual point is 
unity, we can determine the mean to arbitrarily high preci- 
sion if we have a large enough number of points, (Var(/i) — > 
0, in the limit N — > oo). We can therefore differentiate be- 
tween two samples with slightly different means (Sfi <C 1), 
provided we have enough points in each sample. In this ex- 
ample, each point has very little strength (since the difference 
in the means is much smaller that the variance) but in the 
limit of large samples we can detect a small difference in the 
means to high significance. 

The estimation of cosmological parameters from mea- 
surements of the power spectrum is analogous to the above 
example. The CMB satellites will provide, via Ce, measure- 
ments of 1000-3000 nearly statistically independent num- 
bers from which we wish to determine a much smaller num- 
ber (perhaps 10-20) cosmological parameters. Since we have 
a large number of independent points, it is possible to mea- 
sure parameter differences to high significance even though 
the strength of each individual point is small. Furthermore, 
since the number of cosmological parameters that we wish 
to measure is much smaller than the number of independent 
multipoles, we must measure each Ce to an accuracy that 
is much smaller than the cosmic variance; correlated sys- 
tematic errors between the measured Ce coefficients which 
introduces a difference in \ of the order the number of pa- 
rameters can lead to a significant bias in the cosmological 
parameter estimates, the exact bias depending on the na- 
ture of the correlations. We will expand on these points in 
the rest of this Section. 

Figure [l^ is designed to illustrated the concepts of 
strength and significance as applied to measurements of Ce- 
The dotted curves in the upper panels of the figure show Ce 
for a spatially flat CDM model with the parameters listed in 
the Figure heading. The solid lines show the power spectrum 



of the spatially flat target model defined in the previous Sec- 
tion. The parameters for the three examples in this figure 
were chosen so that the Ce curves would be distinguishable 
by OMAP at about the 2a level; in other words, we have 
chosen three pairs of spatially flat CDM models that would 
be difficult to distinguish with OMAP alone even if there 
were no significant systematic errors in the measurements. 
The three examples have been chosen so that: in Figure 
|l^(a) the baryon and CDM densities each differ by about 
12% from those of the target model, while the inflationary 
parameters n 3 , nt and r are almost the same as those of the 
target model; in Figure |l2|(b) , the baryon density differs by 
24% from that of the target model while uj c differs by 6%; 
in Figure [jj|(c) the CDM density differs by 18% from that 
of the tar get model while ut, differs by 4%. In the models 
of Figures |l2](b) and [l2|(c), the large changes in the densi- 
ties have been balanced by relatively large changes in the 
inflationary parameters H. The points in each figure show a 
simulated Ce spectrum for each experiment, assuming the 
dotted curve for Ce, and with errors computed from equa- 
tion (pa|). The lower panels show the residuals relative to the 
target model. 

The number labelled \ 2 m each panel gives the like- 
lihood ratio x 2 ~ —2ln(jC/jCmax), computed from equa- 
tion (j|). (In computing \ 2 w e used i ma x = 1000 for OMAP, 
£ma X = 1500 for CMAP and £ max = 2800 for Planck. How- 
ever, the value of x 2 is insensitive to l ma x since the errors 
in Ce increase exponentially at high multipoles.) The spa- 
tially flat models that we have considered are defined by 7 
parameters, thus a 95% rejection of a model requires a \ 2 
of m 14. This is similar, by construction, to the values of 
X 2 listed for OMAP. For CMAP, and particularly, Planck, 
the values of \ 2 are much larger, demonstrating their su- 
perior ability to distinguish between cosmological models. 
An experiment such as Planck can distinguish between two 
Ce curves at a high significance level even if the differences 
between them are much smaller than the cosmic variance; 
in fact the power spectra and residuals for a pair of mod- 



§ Note that to test inflationary models at interesting levels, we 
would certainly want to determine the spectral index n s to better 
than ±0.05 (see e.g. Liddle 1997 and references therein) 
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(a) cJ b = 0.Olll, cj o = 0.209, h = 0.541, n B =0.995, n t =0.036, r=0.425 
OMAP CMAP PLANCK 




100 300 500 1000 2000 200 500 1000 2000 200 500 1000 2000 
III 

(b) t^ b = 0.0155, cj c = 0.251, h = 0.538, n s =1.0B3, n t =-0.334, r=0.619 
OMAP CMAP PLANCK 




100 200 500 1000 2000 200 500 1000 2000 200 500 1000 2000 

III 
(c) t^ b = 0.0130, cj c = 0.194, h = 0.608, n s =1.052, n t =-0.160, r=0.633 
OMAP CMAP PLANCK 




500 10002000 
I 



500 1000 2000 
I 



500 1000 2000 
I 



Figure 12. Temperature power spectra for three spatially flat models that are distinguishable by OMAP from the target model at about 
the 2a level. The parameters for each model are listed at the top of each figure and the temperture power spectrum computed for these 
parameters is shown by the dashed lines in the upper panels. The power spectrum of the target model is shown by the solid lines. The 
filled circles and error bars show a random realization of Ce for each model and la errors assuming the experimental characteristics of 
OMAP, CMAP and Planck as listed in Table ^. The points in the lower panels show the residuals with respect to the target model. 
The number in each panel gives \ 2 P er degree of freedom for each experiment. The points have been averaged in 100 equally spaced 
logarithmic intervals in I over the range 3 < I < 3000. 
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els with x 2 ~ 14 would be indistinguishable plotted on the 
scales of Figure [T2J 

Small, correlated errors in Ce can thus introduce biases 
in estimates of cosmological parameters. The degree of bias 
depends on the size of the errors and the extent to which the 
correlated component resembles the shape of the derivatives 
of Ce with respect to cosmological parameters. The analysis 
of correlated errors in Ce can therefore be viewed as an- 
other aspect of the analysis of degeneracies in cosmological 
parameter estimation. There are many possible sources of 
correlated errors in Ct. The most obvious are errors caused 
by inaccurate subtraction of emission from the Galaxy and 
extragalactic point sources. However, there are many other 
more instrument specific sources of error, for example, in a 
'total-power' experiment like Planck, the scanning strategy 
together with instrumental '1//' noise can introduce sys- 
tematic errors in maps of the CMB (sometime called 'strip- 
ing', Wright 1996, Bersanelli et al. 1996, Tegmark 1997a, 
Delabrouille 1998). The resulting effect of scanning on Ce 
depends on the details of the experiment and map-making 
methods, but for a fairly wide class of possibilities leads to a 
contribution to Ce that varies as 1/1 (Tegmark 1997a). The 
scanning errors in the CMB maps should be much smaller 
in a differential experiment, such as MAP. Fortunately, fore- 
ground and scanning errors are unlikely to produce wiggles 
in Ce that can distort or mimic the Doppler peak struc- 
ture at high multipoles. Instead, we would expect smoothly 
varying errors to correlate with the low-order principal com- 
ponents. 

To illustrate the effects of correlated errors we have con- 
structed the following simple models 

£(£+l)AC[ 



2tt 

and 

£(£ + l)ACf 
2tt 



= A £(£ + 1) 



1000 

(5 + ey 



+ B 



i *(*+!) 
1000 2 



= K 



a £ K £ Z {£+1) 
(£ 2 + £ 2 K ) 3 / 2 



(32) 



(33) 



The first model (which we denote 'F') is designed to mimic 
the power spectrum of contaminating foregrounds. The first 
term is a rough approximation to the power spectrum of the 
major sources of Galactic emission, free-free, synchrotron 
and dust emission (Low and Cutri 1994; Guarini, Melchiorri 
and Melchiorri 1995, Tegmark and Efstathiou 1996, Bouchet 
et al. 1997) . The second term represents the white noise 
from unresolved extragalactic point sources. The analysis of 
Bouchet et al. (1997) suggests the following values for A and 
B (in llK) at 90GHz and 217GHz: 



A = 0.36, 
A = 1.5, 



B = 1.07, 
B = 1.39, 



(90GHz), 
(217GHz). 



(34a) 
(34b) 



The values of A in (34a 34b) are for the cleanest regions of 
the sky. The second model (which we denote 'S') is based on 
the analysis of Tegmark (1997a) and is designed to approxi- 
mate correlated scanning errors. Values for the constants K 
and Ik are extremely sensitive to the experimental details 
but plausible values for a 'total-power' experiment may be 
£ K « 20 with K in the range 0.004 — O.O^K (see Tegmark 
1997a, figure [y]). The models of equations ([32] ) and ( |33| ) are 
illustrated in Figure Jl3J. 

Figures ^ and [15| shows the biases in a number of cos- 
mological parameters for various choices of the parameters 



a. 



is 



100 



10 



0.1 



- v ^ F (217GHz) 
F (90GHz) y ' ' 



S (K=0.04/iK) 



10 



100 



1000 



Figure 13. The solid line shows the CMB power spectrum of the 
spatially flat target model plotted with a logarithmic ordinate. 
The two dashed lines show the power spectra of the foreground 
model of equation ( fj^ ) with the coefficients A and B set to the val- 
ues appropriate for foregrounds at 90 and 217 GHz as estimated 
by Bouchet et al. (1927). The dotted line shows the scanning error 
model of equation fej) with £pc = 20 and K = 0.04fiK. 



A, B and K. To compute these figures, we have generated 
Monte-Carlo realizations of the CMB power spectrum, to 
which we have added the systematic error models of equa- 
tions ( pS2| ) and (pS3|). We then compute values of the cosmo- 
logical parameters by fitting the first order expansion to Ce 
(equation p^[ ) as in the computations of Figures 0a and Qd. 
The bias in each parameter has been divided by the expected 
standard deviation in each parameter to illustrate when sys- 
tematic errors become larger than the random errors. The 
results are encouraging. For OMAP and CMAP, the tensor 
parameters are significantly affected for values of A ^ 1/iK 
and the remaining parameters for values A £ 3(iK. OMAP 
and CMAP are unaffected by point source contributions 
with B 10/iK because they are insensitive to the power 
spectrum at high multipoles where point sources dominate. 
For Planck, the systematic biases become significant for val- 
ues of A ^ 1/iK and for B Ji, 0.5^iK. Our analysis therefore 
suggests that at the highest frequencies, and in clean regions 
of the sky, MAP should provide unbiased estimates of cosmo- 
logical parameters without any correction for foregrounds. 
For Planck, removal of foregrounds is likely to be required to 
achieve unbiased results. However, with the wide frequency 
coverage and high sensitivity of Planck it should be possi- 
ble to remove foreground contributions to Ce to an accuracy 
that is well below the cosmic variance errors (Tegmark and 
Efstathiou 1996). Even without any frequency information, 
it would be possible to test for biases by including paramet- 
ric models of the systematic errors in the likelihood analysis. 

Figure ^ shows that for Planck, low-level scanning er- 
rors could potentially lead to biases in cosmological parame- 
ters. These biases could probably be removed adequately by 
including a parametric model of the scanning errors derived 
from simulations such as those described by Delabrouille 
(1998). The example of Figure does, however, show that 
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0.1 



CMAP 



PLANCK 




10 



A OiK] 





B [fiK] 

Figure 14. The absolute bias of various cosmological parameters 
divided by their standard errors plotted against the parameter 
A representing Galactic foregrounds, and B representing point 
sources (equation Eil) . 



scanning errors may be a serious concern for the analysis of 
Planck data. 



6 EXTERNAL CONSTRAINTS AND 
PARAMETER DEGENERACIES 

In practice, observations of the CMB anisotropies will be 
supplemented by other astronomical observations which can 
be used to remove some of the parameter degeneracies. Some 
of these are described in this Section. Indeed, such con- 
straints have already been applied in some analysis of cos- 
mological parameters. For example, Webster et al. (1998) 
describe joint likelihood analysis of CMB measurements 
and observations of large-scale structure. White (1998) and 
Tegmark et al. (1998) describe how CMB measurements 
can be combined with distances derived from Type la su- 
pernovae {e.g. Perlmutter et al. 1997, 1998) to remove the 
geometrical degeneracy. 

6.1 Constraints on the Hubble Constant 

If the matter content parameters u m and ujd are well con- 
strained by CMB observations, a constraint on the Hubble 
constant converts directly into a constraint on Qa- For our 
spatially flat target model, the constraint is given by equa- 
tion (n2) and generally, 



Figure 15. The absolute bias of various cosmological parameters 
divided by their standard errors plotted against the parameter K 
representing scanning errors according to equation p3|. 



, 5Ho 
' Ho ■ 



(35) 



A determination of Oa to an accuracy of ~ 0.1 should be 
possible if the systematic and random errors on the Hubble 
constant can be reduced to ^ 5%, as seems feasible (see e.g. 
Freedman et al. 1998). 



6.2 Constraints on the age of the Universe 

The age of the Universe is give by the integral, 



t = 9.8 Gy 



a da 



[u> m a + Uka 2 + LOAa 4 ] 1 / 2 



(36) 



(see e.g. Peebles, 1993, equ 13.9). If the matter content uj m is 
well constrained by the CMB anisotropies, then limits on the 
age of the Universe can be used to constrain the geometrical 
degeneracy in the uia - u>k plane. This is illustrated in Figure 
|is| . However, since the slope of the geometrical lines are so 
steep, the limits on uja are sensitive to the parameters of 
the target model. In the example shown in Figure [if], an age 
constraint of to ^ 14 Gyr (e.g. Chaboyer et al. 1996) would 
lead to a relatively weak limit of uja 0.7 in a spatially fiat 
universe (1Z = 2), but to tighter limits in a universe with a 
larger value of cuk- 



6.3 Constraints on Large-Scale Clustering 

Measurements of galaxy clustering can provide constraints 
on the shape and amplitude of the galaxy power spectrum. 
These are often characterised by the rms amplitude of the 
galaxy fluctuations in spheres of radius 8/i _1 Mpc, (as) g , and 
a shape parameter F (see Efstathiou, Bond and White 1992). 
When combined with CMB anisotropy measurements, ob- 
servation of galaxy clustering can constrain the distribu- 
tion of mass fluctuations relative to galaxies. This informa- 
tion, together with constraints on redshift space anisotropies 
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Figure 16. Analogue of Figure |ljb for a model with u> m = 0.075. 
The solid line show constant values of Ti in the u\— lu^ plane and 
the dashed lines show lines of contours of equal age listed in 10 
Gyr units. 



and cosmic velocity fields, can be used to improve the es- 
timates of many of the cosmological parameters described 
here and including the epoch of reionization (see for example 
Tegmark 1997b). In the future, observations of weak gravita- 
tional lensing may provide more direct estimates of the mat- 
ter power spectrum (Blandford et al. 1991, Miralda-Escude 
1991, Kaiser 1992, see Seljak 1997b for a recent discussion). 
If we restrict to CDM models with low baryon fractions, then 
the power spectrum shape parameter F is related to the mat- 
ter density according to F ~ Q. m h = u) m /h (see Eisenstein 
and Hu, 1998, equations 30 and 31 for a fitting function for 
r that includes the baryon density and associated scale de- 
pendence). A measurement of V from large-scale structure 
measurements can thus be combined with an accurate de- 
termination of u) m from the CMB anisotropies to yield an 
estimate of the Hubble constant, and hence to constrain the 



cosmological constant as described in Section 3.1 



6.4 Constraints on the Geometry from Type la 
Supernovae 

The distances of Type la supernovae at redshifts z 0.5 are 
already providing strong constraints on the geometry of the 
Universe {e.g. Perlmutter et al. 1997, 1998). Such measure- 
ments provide a potentially powerful method of removing 
the geometrical degene racy in the CMB measurements. This 
is illustrated by Figure |l6], which shows the CMB geometri- 
cal degeneracy lines in the Ha — fife plane plotted as in Fig- 
ure |l|a. The dashed lines show constant values of the lumi- 
nosity distance at z = 0.5^], approximating the slope of the 
error ellipses in this plane derived from supernovae measure- 
ments. To first order, standard candles such as Type la su- 
pernovae yield near degenerate constraints along the dashed 

" This is approximately the median redshift of the distant su- 
pernovae samples, (Ellis, private communication). 



Figure 17. The solid lines show the degeneracy lines of con- 
stant 1Z as plotted in Figure Ek. The dashed lines show lines of 
constant luminosity distance at a redshift of z = 0.5 in steps of 
50h~ 1 Wlpc, extending from 1700/i _1 Mpc at the bottom of the 
figure to 2100h _1 Mpc at the top. 



lines. As has been noted as well by White (1998), Tegmark 
et al. (1998) and others, these intersect almost orthogonally 
with the CMB geometrical degeneracy lines. Type la super- 
novae thus offer an extremely powerful way of breaking the 
geometrical degeneracy, provided systematic errors (which 
move the supernova maximum likelihood almost along the 
geometrical degeneracy lines) can be shown to be small. 



7 CONCLUSIONS 

High precision observations of CMB anisotropies promise a 
revolution in our knowledge of fundamental cosmological pa- 
rameters. In this paper we have analysed the degeneracies 
between cosmological parameters determined from CMB ex- 
periments, using the MAP and Planck satellites as an indi- 
cation of what might be achieved in the next decade. Our 
conclusions are as follows: 

[1] The geometrical degeneracy between £Ia and fife for mod- 
els with identical matter content, uib, u) c , ■ ■ - is nearly exact 
and cannot be broken from observations of the linear CMB 
anisotropies alone. Breaking of this degeneracy requires ad- 
ditional data, such as distances to Type la supenovae, or 
the detection of small non-linear effects such as gravitational 
lensing of the CMB. 

[2] Small numerical errors in calculating the CMB power 
spectrum Ci can accumulate in computations of the Fisher 
matrix leading to spurious breaking of near exact degenera- 
cies. These numerical errors can be avoided by rotating to 
new va riab les defined by the degeneracies as described in 
Section 



3.3 



[3] The Fisher matrix for a particular experiment defines a 
set of principal components, i. e. orthogonal linear combina- 
tions of the cosmological parameters Si defining the theoret- 
ical CMB power spectrum Ci . Restricting to a spatially flat 
universe and an idealised set of 7 cosmological parameters, 



© 0000 RAS, MNRAS 000, 000-000 



24 G. Efstathiou and J.R. Bond 



we find that most of the variance in estimates of the Si from 
the MAP satellite comes from the two lowest order (i.e. least 
well determined) principal components. 
[4] The degeneracies between poorly determined and well 
determined parameters limits the accuracy of the Fisher 
matrix. An exact likelihood analysis for a MAP-type ex- 
p erime nt shows that the Fisher matrix can overestimate the 
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as cub by factors of two or more. 

[5] For a Planck-type experiment, the principal components 
are almost identical to physical variables. The five highest 
order principal components couple strongly to u>b, (±>c, Q, n s 
and wa respectively. The two lowest order components cou- 
ple strongly to the parameters r and nt defining the tensor 
component. 

[6] For a MAP-type experiment, the errors on the matter 
densities Ub and uj c are sensitive to theoretical constraints 
on the tensor and scalar fluctuation spectra (r, nt and n s ) 
and hence on the nature of the inflationary model. The vari- 
ances on uib and uj c from a Planck-type experiment are in- 
sensitive to the details of inflation, provided the tensor and 
scalar fluctuation spectra are characterised by weakly scale- 



dependiuiL puwei-laws. 

[7] The lowest order principal components for an experiment 

sampling low multipoles (£ 300 in a spatially flat universe) 
can be understood in terms of the position and height of the 
first Doppler peak. 

[8] Correlated errors in Ce that are much smaller than the 
errors from cosmic variance can introduce significant biases 
in the estimates cosmological parameters. We have anal- 
ysed simple models for correlated errors in Ce arising from 
Galactic foregrounds and extragalactic point source and 
from scanning errors. If our simple models are correct, fore- 
grounds should not introduce large biases in cosmological 
parameters estimated from CMB measurements at frequen- 
cies in the range ~ 100-200GHz. However, there may be 



foregro md components that are not included in our model 
and that dominate in this frequency range (e.g. Leitch et 
al. 1997). In this case, subtraction of the frequency depen- 
dent foregrounds from the primary anisotropies to high ac- 
curacy may be necessary to derive cosmological parameter 
estimates from the CMB. Scanning errors are another po- 
tential source of biases in cosmological parameter estimates. 
Our models suggest that scanning errors may be problematic 
for Planck and a more detailed analysis is evidently required 
before the scanning strategy is frozen. 

[9] Observations of CMB anisotropies can be supplemented 
with more conventional astronomical measurements, such as 
the large-scale mass distribution, estimates of the age of the 
Universe and distances to Type la supernovae, to break the 
geometrical degeneracy. 
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